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Abstract. A 12 C 32 S, 13 C 32 S, 12 C 34 S, and 12 C 33 S J =2-1 
line survey has been made to study interstellar 32 S/ 34 S 
and 34 S/ 33 S ratios from the galactic disk. The four CS 
isotopomers were detected in 20 star forming regions with 
galactocentric distances between 3 and 9 kpc. From a com- 
parison of line velocities, the C 33 S J=2-l rest frequency is 
~ 250 kHz below the value given in the Lovas (1992) cat- 
J>^alog. Taking 12 C/ 13 C ratios from Wilson & Rood (1994) 
i Ctf and assuming equal 12 C 32 S and 13 C 32 S excitation temper- 
atures and beam filling factors, 12 C 32 S opacities are in the 
range 3 to 15; average 32 S/ 34 S and 34 S/ 33 S isotope ratios 
are 24.4 ± 5.0 and 6.27 ± 1.01, respectively. While no 
systematic variation in the 34 S/ 33 S isotope ratio is found, 
the 32 S/ 34 S ratio increases with galactocentric distance 
when accounting for the 12 C/ 13 C gradient of the galactic 
disk. A fit to the unweighted data yields 32 S/ 34 S= 3.3 ± 
0.5(Z>Gc/kpc) + 4.1 ± 3.1 with a correlation coefficient 
of 0.84. Since the interstellar sulfur (S) isotopes are syn- 
thesized by oxygen burning in massive stars, consequences 
— for nucleosynthesis and models of chemical evolution are 
briefly discussed. 
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1. Introduction 

All elements heavier than hydrogen are believed to be 
synthesized by nuclear processes. The standard Big Bang 
model can explain the bulk of the elemental abundances 
up to 7 Li, while heavier nuclei must be synthesized in stars 
(e.g., Wilson & Rood 1994). Thus accurate measurements 
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of elementary and isotopic abundances are needed to gain 
insight into the physical processes in the interior of stars. 

While information on stellar and interstellar abun- 
dances can be obtained from optical spectroscopy, such 
studies of the atomic gas phase are mainly confined to the 
solar neighborhood and do not allow to discriminate be- 
tween various isotopic constituents. This can only be done 
by observing molecular isotopomers. Such studies, carried 
out at cm- and mm-wavelengths, have the additional ad- 
vantage that they are probing optically obscured regions 
across the entire Galaxy, thus providing information on 
its 'chemical' history. While the isotopic abundances of C, 
N, and O provide detailed information on CNO and he- 
lium burning (e.g., Wilson & Matteucci 1992; Henkel & 
Mauersberger 1993; Henkel et al. 1994b; Wilson & Rood 
1994; Kahane 1995), sulfur (S) is a useful tool to study 
the late evolutionary stages of massive stars. 

Among the sulfur bearing molecules, CS (carbon 
monosulfide) is the most suitable species: it is ubiquitous 
in the interstellar medium (ISM), its main isotopic lines 
are strong, and its lowest rotational transitions are ob- 
served at frequencies not requiring extremely favorable 
weather conditions. The three most prominent sulfur iso- 
topes, 32 S, 34 S, and 33 S have solar system abundance ra- 
tios of 22.5 : 1 : 0.18 (e.g., Bogey et al. 1981; Kahane 
et al. 1988). These abundance ratios were also found in 
the late-type giant IRC+10216 (Wannier 1980; Kahane et 
al. 1988). The 32 S/ 34 S ratio was already measured with 
a beam size of 130" toward star forming regions of the 
galactic disk and center (Frerking et al. 1980). During the 
last fifteen years, larger telescopes and more sensitive re- 
ceivers have led to a higher quality of observational data. 
Also, interstellar 12 C/ 13 C ratios which are needed to ob- 
tain reliable 32 S/ 34 S values (see Sect. 4) are now much 
better known. Although C 33 S was detected by Sutton et 
al. (1985) and Cummins et al. (1986), no systematic study 
of interstellar 34 S/ 33 S ratios has yet been reported for the 
ISM. We thus observed the J=2-l lines of 12 C 32 S (here- 



C 34 S), and 12 C 33 S (hereafter C 33 S) toward a number of 
star forming regions. The 13 CS, C 34 S, and C 33 S lines allow 
to measure the relevant sulfur isotope column densities, 
while 12 CS and 13 CS are used to estimate the opacity of 
the 12 CS J =2-1 line. 



2. Observations 

The measurements were made in September 1993 and May 
1994 using the SEST 15-m telescope at La Silla, Chile. At 
the frequencies of the J=2-l transitions of 12 CS, 13 CS, 
C 34 S, and C 33 S (97980.986, 92494.299, 96412.982, and 
97172.086 MHz, respectively; Lovas 1992) the beamwidth 
was ~ 53". We employed a Schottky receiver with system 
temperatures of order 500 K on a main beam brightness 
temperature scale (Tmb)- The backend was an acousto- 
optical spectrometer with 2000 contiguous channels and 
a channel separation of 43kHz (0.13 - 0.14kms _1 ). All 
measurements were carried out in a dual beam switching 
mode (switching frequency 6 Hz) with a beam throw of 
11'40" in azimuth. Integration times varied from 4 min- 
utes for 12 CS to 40 minutes for C 33 S. 

The pointing accuracy, obtained from measurements 
of SiO maser sources, was better than 10". Using the soft- 
ware package CLASS for data reduction, all spectra were 
converted to a Tmb scale, adopting a main beam efficiency 
of »7mb=0.74. Within a single observing period, line tem- 
peratures could be reproduced to 10%. Comparing C 34 S 
data from Sept. 1993 and May 1994, the Tmb scales are 
identical within 15%. C 34 S and C 33 S spectra from Orion- 
KL, taken in January 1995, also show only small devi- 
ations (<10%) from the results displayed in Fig. 1 and 
Table 1. This is consistent with the known long term sta- 
bility of the 3 mm Schottky receiver which allowed to re- 
produce line intensities from NGC 4945 over a period of 
several years (cf. Henkel et al. 1990, 1994a; Whiteoak et 
al. 1990; Bergmann et al. 1992; Israel 1992; Dahlem et al. 
1993). Furthermore we obtained in May 25, 1994, C ls O 
and C ir O SEST spectra and in June 1, 1994, C ls O and 
C ir O 30-m IRAM spectra from M 17 SW (ai 950 = 18h 17 m 
30 s , 61950 = -16°12'58"). The SEST C 18 0/C 17 line in- 
tensity ratios obtained in a frequency switching and a dual 
beam switching mode (see above) were 3.43 and 3.50, re- 
spectively. The 30-m IRAM maps consisted of 16 positions 
with 20" spacing centered on the SEST position; these 
data were obtained in a position switching mode with 
the reference at an offset (Aa, A<5)= (20', 0). Convolv- 
ing the IRAM spectra to a full width to half power beam 
size of 48" (the resolution of the SEST data), we obtain 
I(C 18 0)/I(C 17 0) = 3.47. For a given position, the error 
in the individual J(C 18 0)/J(C 17 0) ratios due to noise is 
~ 10% at the 15-m SEST and 30-m IRAM telescopes. The 
good agreement between the SEST and IRAM data lends 
further confidence to the results of this study. 



Wood & Churchwell (1989) have shown that the colors 
of known ultracompact (UC) HII regions, i.e. recently 
formed O and B stars still embedded in their natal molec- 
ular cloud, are narrowly confined in the far infrared (FIR) 
color-color plots. They further showed that there are few 
other objects which share the FIR colors of the UC HII 
regions. Applying the two color criterion, it is then possi- 
ble to identify embedded O and B stars anywhere in the 
galactic disk. From a complete sample of 74 IRAS sources 
with 8 < —30° and 5ioo^m>3000 Jy, we selected 20 ob- 
jects with T MB (C 34 S) > 0.85 K for further 12 CS, 13 CS, and 
C 33 S measurements. Besides these IRAS sources, Orion- 
KL and two positions toward the north-eastern part of the 
prominent star forming region NGC 6334 (radio source F 
(Rodriguez et al. 1982) and FIR-I (Straw & Hyland 1989)) 
were also observed. The distances of the southern IRAS 
objects to the galactic center, -DqCj have been estimated 
with the following model: 

®i£2£l = 1.0074 ( 0382 + 0.00698 . (1) 

©(-Dgc) is the rotational velocity of a source at galacto- 
centric distance -Dqc! ©o = 220kms _1 (Brand 1986) and 
D & = 8.5 kpc is the galactocentric distance of the Sun. 

4. Results 

The measured CS profiles are shown in Fig. 1 after a first 
order baseline subtraction and a removal of the image side- 
band 12 CS feature in the C 33 S spectra. The weak features 
sometimes seen at the blueshifted and redshifted sides of 
the C 33 S profile could not be identified (cf. Lovas 1988). 
From gaussian fits to the line parameters of NGC 6334 A, 
the frequency offsets of these weak features relative to 
C 33 S are +3.4 ± 0.1 and -2.3 ± 0.1 MHz, respectively 
(for z>"(C 33 S), see Sect. 2). The corresponding image side- 
band frequencies are ~ 97991 and 97997 MHz. The CS line 
parameters obtained from gaussian fits are summarized in 
Table 1. 

There is no significant difference in the radial veloci- 
ties between 12 CS, C 34 S, and 13 CS (A «lsr( 13 CS, 12 CS) 
= 0.15 ± 0.32kms- 1 ; A v L sr( 13 CS, C 34 S) = 0.00 ± 
0.21kms _1 ). H owever, the C 33 S velocities are different. 
In all 20 sources, it is C 33 S which has the highest velocity 
(A v L sr(C 33 S, 12 CS) = 0.88 ± 0.27kms _1 ; A v L sr(C 33 S, 
C 34 S) = 0.73 ± 0.29kms _1 ; A v L sr(C 33 S, 13 CS) = 0.73 
± 0.29 km s" 1 ). From our Orion-KL data taken in January 
1995, A v L sr(C 33 S, C 34 S) = 0.46 ± 0.07 kms -1 . This is 
consistent with the line parameters displayed in Table 1. 
We conclude that the C 33 S J=2-l frequency measured by 
us is ~ 250 kHz higher than the calculated frequency given 
in the Lovas (1992) catalog. According to Bogey et al. 
(1981), the central and strongest group of hyperfine com- 
ponents is at v = (97171.840 ± 0.030) MHz. This is 246 



Fig. 1. 12 CS, 13 CS, C 34 S, and C 33 S spectra from our sample of 20 southern sources 



Table 1. Line parameters from unsmoothed J=2-l spectra of the four measured CS isotopomers 
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catalog and agrees well with our observational result. 

The main isotopomer, 12 CS, is known to be moder- 
ately optically thick (e.g., Frerking et al. 1980; Linke & 
Goldsmith 1980). Therefore, the 32 S/ 34 S isotope ratio can- 
not be determined from the observed I( 12 CS)/I(C 34 S) 
ratio. With the assumption that C 34 S and 13 CS are 
optically thin, however, 32 S/ 34 S can be estimated from 
I( 13 CS)/I(C 34 S), if the 12 C/ 13 C ratio is known: 

32g 12 C J(13 C 32 S ) 

34 S ~ 13 C I( 12 C 34 S)' ^ 

Assuming equal excitation temperatures and beam filling 
factors for 12 CS and 13 CS, the 12 CS J =2-1 peak opacity 
r ( 12 CS) can be determined from 

T MB ( 12 CS) l-e-(' 2cs ) 

Tmb( 13 CS) ~ 1-e-r ( 12 cs)/B' K i3 C ' W 

12 C/ 13 C ratios have been obtained from H2CO (e.g., 
Henkel et al. 1982, 1985) and C ls O (Langer & Penzias 
1990) throughout the galactic disk. In the ISM, the carbon 
isotope ratio is thus the best known among the various 
CNO isotope ratios. Wilson & Rood (1994) obtain from a 
fit to these data 

12 C 

— = (7.5 ± 1.9) (£> GC /kpc) + (7.6 ±12.9). (4) 

For Orion-KL, 12 C/ 13 C ratios of 79 ± 7 and 67 ± 3 were 
determined from C ls O (Langer & Penzias 1990). Because 
of possible errors in a single measurement, fractionation 
(see Sect. 5.2), and selfconsistency, Eq. 4 was also applied 
to the Orion-KL data. 

While the 32 S/ 34 S isotope ratio cannot be determined 
from CS alone, the 34 S/ 33 S isotope ratio is directly ob- 
tained from 

^S J(C 34 S) 

33 S ~ I(C 33 S)' (5) 

Eq. 5 requires that C 34 S is optically thin (see Sect. 5.1). 

In Table 2, galactocentric distances, 12 CS J =2-1 peak 
opacities, and 12 C/ 13 C, 32 S/ 34 S, and 34 S/ 33 S isotope ra- 
tios are derived from Eqs. 1-5 and Table 1 for our sample 
of sources. The average 32 S/ 34 S and 34 S/ 33 S ratios are 
with 24.4 ± 5.0 and 6.27 ± 1.01 (the errors are the stan- 
dard deviations of the mean) consistent with the solar 
system values of 22.5 and 5.5. 

5. Discussion 

5.1. Photon trapping and line saturation 

To derive the quantities displayed in Table 2, unsaturated 
13 CS, C 34 S, and C 33 S J=2-l lines as well as equal excita- 
tion temperatures (T ex ) for all isotopomers have been as- 
sumed. In view of significant saturation and photon trap- 
ping in the main isotopic species, the validity of the two 
assumptions has to be examined in more detail. 



temperatures (T ex ) relative to the optically thin limit, can 
be estimated when calculating the statistical equilibrium 
populations of the various CS isotopes in the J=l and 
2 states. We have performed such calculations using the 
LVG (Large Velocity Gradient) approximation for a cloud 
of spherical geometry (e.g., Scoville & Solomon 1974). 
With collisional rates from Green & Chapman (1978), the 
populations for the 13 lowest rotational CS states have 
been determined. Table 3 provides some characteristic re- 
sults for a kinetic temperature Tk; n = 30 K (cf. Padman et 
al. 1985), n(H 2 ) = 10 5 cm -3 (cf. Mauersberger & Henkel 
1989), and isotope ratios of 1 : 23, 1 : 60, and 1 : 126.5 
relative to the main species. Note that saturation in C 34 S 
will decrease the I(C 34 S)/I( 13 CS) and I(C 34 S)/I(C 33 S) 
line intensity ratios below the corresponding abundance 
ratios. This effect is however significantly reduced by pho- 
ton trapping. As exemplified by Table 3, the line intensity 
ratios should provide reliable isotope ratios up to r (C 34 S) 
~ 1, where the error due to saturation becomes ~ 20% or 
less. In view of the r max values given in Col. 5 of Table 2, 
this requirement appears to be fulfilled for most sources. 

A direct test of the significance of C 34 S line saturation 
to our data can be made by plotting the J( 13 CS)/J(C 34 S) 
and J(C 33 S)/J(C 34 S) line intensity ratios as a func- 
tion of I(C 34 S)/I( 12 CS). Sources with large values of 
I(C 34 S)/I( 12 CS) are more likely to be optically thick in 
C 34 S, which should also enhance the I( 13 CS)/I(C 34 S) and 
I(C 33 S)/I(C 34 S) ratios. As can be seen from Figs. 2a and 
b, no correlation is evident. This indicates that for most of 
the disk sources C 34 S saturation is not significant, while 
a few optically thick sources cannot be excluded. 

Another way to estimate the degree of saturation in 
the C 34 S J=2-l lines is a comparison of linewidths. We 
obtain for all 20 sources A (Av 1/2 ( 12 CS,C 34 S)) = 1.05 ± 
0.74kms _1 , A (Av 1/2 (C 34 S, 13 CS)) = 0.19 ± 0.41 kms" 1 , 
and A (Av 1/2 (C 34 S,C 33 S)) = -0.04 ± 1.03kms _1 . For 
the five spectra with the highest quality C 33 S pro- 
files (toward IRAS 15491, IRAS 16172, IRAS 17059, and 
NGC6334A and B), A (Av 1/2 (C 34 S,C 33 S)) = 0.16 ± 
0.41 kms -1 . We conclude that while the saturation of 
12 CS is indicated by a higher linewidth (Av!/ 2 ( 12 CS) > 
Av 1/2 (C 34 S, 13 CS, C 33 S) in 80% of our sources), there 
is no significant difference between the other isotopomers 
(Av 1/2 (C 34 S) >Av 1/2 ( 13 CS, C 33 S) in only 40% of the 
sources). This result also suggests that C 34 S line satu- 
ration is not important here. 

5.2. Fractionation 

Observed isotope ratios cannot be used for an analysis of 
stellar nucleosynthesis and galactic 'chemical' evolution if 
the results are significantly modified by fractionation. The 
12 C/ 13 C isotope ratio is likely affected by the reaction 

13 C + + 12 C q > 12 C + + 13 C Q + A £ 35K ( 6 ) 



Table 2. Sulfur isotope ratios and 12 CS optical depths 1 



Source 


£>GC b) 

(kpc) 


12 C/ 13 C c) 


Tmb( 12 CS)/T MB ( 13 CS) d ' 


r max ( 12 CS) e ) 


J( 13 


CS)/7(C 34 S) 


32 s/ 34 s f) 


34 S/ 33 S s) 


Orion-KL 


9.0 h ' 


75±21 


18 


69±1 


73 


4.1 





467±0 


015 


35±10 


7 


53±0.45 


IRAS 12320 


7 


j) 


60±19 


10 


53±1 


55 


6.0 





342±0 


019 


21±7 


5 


66±0.32 


IRAS 12326 


7 


1 j) 


61±19 


11 


67±1 


08 


5.4 





433±0 


014 


26±8 


4 


38±0.16 


IRAS 15290 


5 


j) 


45±16 


8 


74±1 


37 


5.4 





438±0 


022 


20±7 


7 


72±0.62 


IRAS 15408 


6 


5 j > 


56±18 


15 


87±2 


36 


3.5 





416±0 


019 


23±8 


8 


87±0.90 


IRAS 15491 


6 


1 j) 


53±17 


7 


99±0 


62 


7.1 





475±0 


010 


25±8 


6 


83±0.25 


IRAS 15520 


6 


2 j) 


54±17 


6 


23±0 


28 


9.4 





536±0 


008 


29±9 


6 


42±0.32 


IRAS 15567 


4 


3 j) 


40±15 


13 


46±2 


03 


2.9 





435±0 


021 


17±7 


6 


08±0.37 


IRAS 16065 


5 


3 j > 


47±16 


6 


87±0 


61 


7.4 





523±0 


014 


25±8 


5 


73±0.26 


IRAS 16164 


5 


4 j) 


48±16 


8 


41±0 


96 


6.1 





507±0 


014 


24±8 


4 


88±0.21 


IRAS 16172 


5 


6 j) 


50±17 


5 


94±0 


29 


9.2 





426±0 


007 


21±7 


5 


46±0.19 


IRAS 16506 


6 


2 j) 


54±17 


9 


52±1 


18 


6.0 





410±0 


021 


22±7 


6 


28±0.53 


IRAS 16562 


7 


j > 


60±19 


9 


63±0 


89 


6.6 





457±0 


012 


27±9 


6 


38±0.38 


IRAS 17009 


6 


5 j) 


56±18 


11 


03±1 


41 


5.3 





357±0 


015 


20±6 


6 


41±0.39 


IRAS 17059 


6 


1 j) 


53±17 


7 


33±0 


47 


7.8 





454±0 


013 


24±8 


5 


24±0.14 


IRAS 17160 


2 


9 j) 


29±14 


8 


34±1 


06 


3.6 





485±0 


019 


14±7 


6 


35±0.36 


NGC6334 A 


7 





60±19 


6 


77±0 


17 


9.6 





500±0 


005 


30±10 


5 


72±0.08 


NGC6334B 


7 


i] 


60±19 


5 


41±0 


20 


12.3 





468±0 


006 


28±9 


6 


68±0.13 


IRAS 17233 


7 


8 j) 


66±20 


4 


89±0 


34 


15.1 





490±0 


009 


32±10 


6 


65±0.18 


IRAS 17257 


6 


3 j > 


55±18 


10 


00±0 


72 


5.8 





461±0 


009 


25±8 


6 


16±0.14 



a) The standard deviations for the Tmb , integrated line intensity, and isotope ratios (Cols. 4-8) do not include systematic cali- 
bration errors which are difficult to assess (for details, see Sect. 2). 

b) The distance of the Sun to the galactic center, D@, is assumed to be 8.5 kpc. 

c) The 12 C/ 13 C ratio is estimated from the distance to the galactic center (see Eq. 4). The error is a = (1.9 2 JJ>gc 2 + 12.9 2 ) 1/2 . 

d) Peak temperature ratios; the uncertainties are derived from the 1<t noise level of the unsmoothed spectra (see Table 1). Note 
that the lines extend over several contiguous channels. The exceptionally large ratio toward Orion-KL is caused by the high 
velocity outflow. 

e) Peak optial depths obtained from Eq. 3. No error is given, since more rotational CS transitions are needed for a successful 
simulation of the physical conditions (see Sect. 5.1). In the case of Orion-KL, the high velocity outflow seen in 12 CS reduces the 
calculated optical depth substantially. 

f) Calculated from Eqs. 2 and 4 

g) Calculated with Eq. 5 

h) see Langer & Penzias (1990) 

i) See Gardner & Whiteoak (1979) 
j) see Eq. 1 



(Watson et al. 1976), which enhances 13 CO in the more 
diffuse C + rich parts of molecular clouds. Eq. 4 is derived 
from both C ls O and H 2 CO data, which are believed to 
bracket the actual carbon isotope ratio (e.g., Langer et al. 
1984; Langer & Penzias 1990). For CS, the corresponding 
reaction 

13 C++ 12 CS > 12 C + + 13 CS+ A £ 25K) (J) 

may be less important. The bulk of the CS emission arises 
from the densest parts of the molecular clouds only (e.g., 
Snell et al. 1984; Mauersberger & Henkel 1989), where 
the relative abundance of C + may be smaller than in the 
more diffuse ISM traced by CO. Also, the kinetic tem- 
peratures to reach a given level of fractionation must be 



smaller for CS than for CO. Most important, however, CO 
can be fractionated and CO is more abundant than CS. 
The CO isotope exchange reaction (Eq. 6) thus enhances 
13 CO and, as it destroys 13 C+, it increases the 12 C+/ 13 C+ 
ratio above the local 12 C/ 13 C ratio. The effective enhance- 
ment of 12 C + over 13 C + implies, that molecules forming 
from C + for which no isotope exchange is possible (e.g. 
CH, H 2 CO, HCN, C 2 H), will show larger 12 C/ 13 C ratios 
than CO and the local ISM. Since CS can be fractionated 
(Eq. 7), it is an intermediate case and we find 

12 CO/ 13 CO < 12 CS/ 13 CS < H^ 2 CO/H^ 3 CO (8) 

for the abundance ratios of CO, CS, and H2CO. In view 
of the moderate differences in the 12 C/ 13 C ratios derived 



Table 3. Computed CS J =2-1 line parameters for Ty ln = 30 K 



ra(CS)/(d*>/dr) a) 


n(H 2 ) 






T b ) 


{ —2 /i — 1\ 

(cm /kms ) 


(cm ) 


( K ) 




( K ) 


3.08 10 14 


10 5 


18.1 


4.52 


14.70 


1.34 10 13 


10 


8.6 


0.62 


2.53 


5.13 10 12 


10 5 


7.7 


0.28 


1.10 


2.43 10 12 


10 5 


7.3 


0.16 


0.55 



a) The ratios of the chosen CS column densities per kms 1 
are 60 : 2.6 : 1 : 0.47 (for 12 CS, C 34 S, 13 CS, and C 33 S, respec- 
tively). 

b) At Tki n = 30 K, the correlation between Tmb values and op- 
tical depths tends to be consistent with Eq. 3 for densities well 
in excess of 10 4 cm -3 . Differences between results obtained in 
Eq. 3 and the LVG code are caused by variations in T ex and by 
the use of Rayleigh- Jeans line temperatures in Eq. 3. The de- 
tailed radiative transfer (LVG) calculations allow to estimate 
excitation temperatures and thus provide those Planck tem- 
peratures which are actually measured. 



from H 2 CO and CO (e.g., Henkel et al. 1982; Langer & 
Penzias 1990), we conclude that the 12 CS/ 13 CS abun- 
dance ratio should be close to the actual 12 C/ 13 C ratio. 

For isotope exchange reactions similar to Eqs. 6 and 7, 
but involving S + , very low temperatures are required. The 
energies of the ground vibrational states differ between 
12 CS and C 34 S (C 33 S and C 34 S) by 7.5 (4) K. This is not 
sufficient to allow for significant fractionation near sites of 
massive star formation. 



5.3. The 32 S/ 34 S gradient 

Fig. 3 displays the 32 S/ 34 S and 34 S/ 33 S isotope ratios as a 
function of galactocentric distance. According to Fig. 3a, 
there is strong evidence for the presence of a 32 S/ 34 S gradi- 
ent. A linear least squares fit to the unweighted data yields 
32 S/ 34 S= 3.3 ± 0.5(D G c/kpc) + 4.1 ± 3.1 with a corre- 
lation coefficient of 0.84. However, the case is not as clear 
as it appears to be at first sight: Most of the sources are 
located within a small Dqc interval ranging from 5.5 to 
7kpc. The apparent increase in the 32 S/ 34 S ratio with in- 
creasing Dqc is thus based on only half a dozen of sources 
outside that Dqc range and an independent confirmation 
of the gradient is clearly needed. 

The 34 S/ 33 S data do not show any hint for a systematic 
variation with galactocentric radius. We suspect that a 
significant part of the scatter in the ratios (almost a factor 
of 2) is due to the small signal-to-noise ratios in some C 33 S 
spectra and blending with the nearby unidentified features 
(cf. Sect. 3). The image sideband 12 CS profile in the C 33 S 
spectra (not shown in Fig. 1; see Sect. 4) is not a problem. 
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Fig. 2. a) 7( 13 CS)/7(C 34 S) and b) 7(C 33 S)/7(C 34 S) as a func- 
tion of 7(C 34 S)/7( 12 CS). In the case of saturated C 34 S J=2-l 
profiles, a positive correlation between the plotted intensity 
ratios is expected. 



5.4- A comparison with other data 

Frerking et al. (1980) presented 12 CS, 13 CS, and C 34 S 
J=2-l data from 14 massive star forming regions. The 
two sources which are in common with our sample are 
Orion-KL and NGC 6334, where a position between the 
two observed by us was chosen. Their J( 13 CS)/J(C 34 S) ra- 
tio toward Orion-KL, 0.36 ± 0.02, is ~ 25% smaller than 
that given in Table 2; their ratio for NGC 6334 is with 0.52 
± 0.03 consistent with our values. Overall, Frerking et al. 
do not find strong evidence for C 34 S saturation and their 
J( 13 CS)/J(C 34 S) ratios from the galactic disk do not vary 
with galactocentric radius. These results are consistent 
with our data. The 32 S/ 34 S gradient, not apparent in the 
I( 13 CS)/I(C 34 S) line intensity ratios, is only seen when 
accounting for the interstellar 12 C/ 13 C gradient across the 
galactic disk. Fifteen years ago, such a systematic varia- 
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Fig. 3. a) 32 S/ 34 S and b) 34 S/ 33 S isotope ratios as a func- 
tion of galactocentric radius. The solar system values at Z>gc = 
8.5 kpc are also given. The relatively large error bars in the 
32 S/ 34 S data are caused by the uncertainty in the 12 C/ 13 C ra- 
tios (see Table 2). For the least square fit to the 32 S/ 34 S data 
displayed by a solid line, see Sect. 5.3. 



tion in the carbon ratio was not well established and was 
thus neglected by Frerking et al. (1980). The absence of 
a significant J( 13 CS)/J(C 34 S) gradient with galactocen- 
tric radius and the agreement in these line intensity ratios 
(0.42 ± 0.16 for the Frerking et al. and 0.45 ± 0.10 for our 
data; see Table 2) does not only provide some additional 
confirmation for the 32 S/ 34 S gradient proposed by us (see 
Sect. 5.3), but also shows that different angular resolutions 
(130" versus 53") do not qualitatively alter the properties 
of the line intensity ratios. 

As mentioned in Sect. 1, C 33 S was first observed by 
Sutton et al. (1985) and Cummins et al. (1986). For a 
30" region centered toward Orion-KL, Sutton et al. obtain 
from the J=5-4 transition I( 12 CS) : I( 13 CS) : I(C 34 S) 
: I(C 33 S) = 291 : 16 : 33.5 : 11.0. While their J=5-4 



corresponding J=2-l line value, their I(C 34 S)/I(C 33 S) 
ratio is only half as large as our value (see Table 2). 
Whether this discrepancy is caused by calibration uncer- 
tainties or other effects (e.g. 'chemical' inhomogeneities) 
remains to be seen. For the J=2-l peak intensity ra- 
tios tward Orion-KL, the survey by Turner (1989) yields 
?mb( 13 CS) : T MB (C 34 S) ~ 0.3 : 0.74. This ratio is interme- 
diate between the integrated line intensity ratios obtained 
by Frerking et al. (1980) and from our data (see Table 2). 

5.5. Sulfur nucleosynthesis 

All the sulfur isotopes discussed in the present paper 
are believed to be of primary nature, synthesized dur- 
ing hydrostatic or explosive oxygen burning. These pro- 
cesses only take place in massive stars and Type la super- 
novae (SNe), with hydrostatic contributions dominating 
in the first, and solely explosive contributions dominating 
in the latter source (Thielemann & Arnett 1985, Weaver & 
Woosley 1993). The relative production ratios of 32 S, 33 S, 
and 34 S in massive stars appear to be rather independent 
of the initial stellar metallicity (Woosley & Weaver 1982). 
According to Thielemann et al. (1994) the present contri- 
bution of Type la SNe to the 32 S production is ~ 40%. A 
comparison of Figs. 4 and 5 of Timmes et al. (1995) indi- 
cates a 20% contribution. While massive stars appear to 
slightly underproduce 33 S but overproduce 34 S compared 
to 32 S (Timmes et al. 1995; this holds equally well for 
SN Ib/c models which we do not distinguish from SN lis 
in the following; cf. Woosley et al. 1995), the situation is 
just the opposite in Type la SNe (Thielemann et al. 1994). 

These nucleosynthesis models do not favor galactic gra- 
dients of the sulfur isotopic ratios. Their primary nature, 
i.e. the independence of their production ratios from the 
initial stellar metal content, does not favor a relation of 
the sulfur isotopic ratios with the (controversial) galactic 
disk metallicity gradient (e.g., Matteucci & Francois 1989; 
Kaufer et al. 1994). A variation of the ratio of galactic 
Type II to Type la SN rates with galactocentric distance 
may also be problematic: This should lead to a signifi- 
cant variation of the 34 S/ 33 S ratio, which is not supported 
by our data (Fig. 3b), although the large scatter in the 
34 S/ 33 S ratios does not exclude small (<50%) variations. 
A galactic 32 S/ 34 S gradient may thus imply that one of 
the basic assumptions in the generally accepted nucleosyn- 
thesis or chemical evolution models needs to be dropped. 
One possibility is that different kinds of Type la SNe, like 
Chandrasekhar mass and sub-Chandrasekhar mass mod- 
els (cf. Woosley & Weaver 1995) yield different sulfur iso- 
tope production ratios, a question which still needs to be 
explored. If so, the data presented here might constrain 
their spatial distribution in the Milky Way. On the other 
hand, even the modelling of hydrostatic oxygen burning in 
massive stars may be generally oversimplified, as recently 
shown by the multi-dimensional calculations of Bazan & 



sulfur production are still to be investigated. 

In any case, the data presented here may turn out 
to provide an important test for the oxygen burning nu- 
cleosynthesis. This is another reason (cf. Sect. 5.3) why 
more observational support for the 32 S/ 34 S gradient is 
still needed. Sources of the northern hemisphere with well 
know distances from the Sun are suitable targets for a 
future study. Such observations would also allow to inves- 
tigate (and to exclude) a Malmquist bias which might be 
caused by different spatial distributions of the various iso- 
topomers (when analysing CNO isotope ratios, this effect 
was found not to be important; the effect is also not appar- 
ent when comparing CS data obtained with 130" and 53" 
beam sizes; see Sect. 5.4). Dark clouds, where ejecta from 
massive stars have not recently contaminated the ISM, 
may yield different ratios. Interesting information can also 
be obtained from the galactic center region. Accepting the 
gradient as shown in Fig. 3a with the solar system ratio 
well below the local interstellar value, we note that the 
results are similar to those obtained for the 14 N/ 15 N ratio 
(Dahmen et al. 1995). 

6. Conclusions 

Having studied the J =2-1 line profiles of 12 CS, 13 CS, 
C 34 S, and C 33 S in 20 prominent star forming regions, we 
obtain the following main results: 

(1) From the CS data and previously determined 
12 C/ 13 C ratios (Wilson & Rood 1994), we obtain a 32 S/ 34 S 
gradient of 3.3 ± 0.5 (-Dcc/kpc) for the inner galactic disk 
out to the solar circle (3kpc < -Dgc< 9kpc). Since only 
a few of our sources are located outside the range 5.5 < 
-Dgc< 7kpc, this result has to be confirmed by additional 
measurements covering the less well sampled galactocen- 
tric distances and including northern sources to firmly ex- 
clude a possible Malmquist bias. A 32 S/ 34 S gradient is not 
expected in current theoretical schemes explaining sulfur 
nucleosynthesis in massive stars and Type la supernovae. 

(2) Models of stellar nucleosynthesis predict that 33 S 
enrichment is dominated by Type la supernovae, while a 
larger part of the 34 S synthesis occurs in massive stars. 
In view of this difference, the absence of a strong galactic 
disk 34 S/ 33 S gradient is an interesting result. 

(3) Assuming equal excitation temperatures and beam 
filling factors, peak 12 CS J=2-l opacities fall into the 
range 3-15. 

(4) The calculated J=2-l C 33 S frequency given in the 
Lovas (1992) catalog appears to be too high. A correction 
of ~ 250 kHz is suggested by the data. 

(5) Additional measurements not only including north- 
ern sources, but also the galactic center region and dark 
clouds may add crucial information on sulfur isotope ra- 
tios and on oxygen burning in Type la supernovae and 
evolved massive stars. 
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